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TECHNICAL NOTE k2k2. 


GENERAL SOLUTIONS FOR FLOW PAST SLENDER CAMBERED WINGS 
WITH SWEPT TRAILING EDGES AND CALCULATION OF 
ADDITIONAL LOADING DUE TO CONTROL SURFACES 
By E. B. Klunker and Keith C. Harder 

SUMMARY 


The slender -wing- type analysis is used to obtain general expres- 
sions for the surface pressure, lift, and rolling moment for cambered 
wings with swept trailing edges. These results are specialized to give 
the additional loading due to the deflection of trailing-edge control 
surfaces and the loading due to a particular type of wing twist. In 
contrast to flat-plate wings, cambered surfaces (or wings with control 
surfaces) support load downstream of the maximum span of the wing where 
the slope of the leading edge is equal to or less than zero. The pres- 
sure distribution and the streamwise lift and rolling-moment distribu- 
tions are discontinuous across the plane normal to the stream direction 
passing through the wing maximum span unless the slope of the wing leading 
edge is continuous and zero at this chordwlse station. A numerical 
example for an aileron control Is included. 


INTRODUCTION 


The flow field about configurations which are slender in the stream- 
wise direction has been shown to assume a two-dimensional character in 
planes normal to the stream velocity. A basically three-dimensional flow 
is thereby approximated by solutions for a two-dimensional crossflow, and 
as a consequence, relatively complex configurations can be analyzed within 
this slender-airplane approximation. Moreover, for configurations with- 
out thickness, the solutions give the first-order approximation for the 
flow at transonic speeds, and the restrictions on the slenderness are not 
so severe as at subsonic or supersonic speeds. The slenderness concept 
has been applied to the analysis of a variety of configurations since the 
Investigation of Jones (ref. l) on the lift of low-aspect-ratio wings. 

The results of many of these studies are presented in references 2 and 5 
together with a discussion of the underlying basis of the theory. Most 
of these Investigations have been restricted to wings with unswept trailing 
edges and therefore have not had to account for the shed vortex sheet in 
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the evaluation of the forces and surface pressures. Some exceptions are 
the studies of Lomax and Heaslet (ref. 4), of Mirels (ref. 5 ), and of 
Mangier (ref. 6). The indirect problem of determining the trailing-edge 
shape of flat-plate wings for a prescribed span loading is treated in 
reference 4, and the lift, pitch, and roll solutions applicable to flat- 
plate wings with swept trailing edges are treated in references 5 and 6. 

The present study is concerned with the development of general solu- 
tions for the flow past arbitrarily cambered slender wings with swept 
trailing edges and the subsequent application of these results to the 
calculation of the additional load due to the deflection of a trailing- 
edge control surface. The engineering usefulness of the solutions for 
arbitrary camber distributions is somewhat limited since the effort 
required to obtain numerical results is not justified in general for 
the slender-wing approximation; their primary value is to show some of 
the general features of the load distribution. The control-surface solu- 
tions as well as solutions for the loading due to pure twist, however, 
can be evaluated without serious difficulty. 

The method of cons time ting the solutions is presented in section 1 
and the general solutions for symmetric and antisymmetric spanwise loading 
are developed in sections 2 and 3* Section 4 presents the additional load 
due to the deflection of a flap or aileron as determined from these gen- 
eral equations. One example is presented to illustrate the method of 
solution for the effect of pure twist on the loading, and calculations 
are presented for the additional loading due to an aileron. A list of 
the more important symbols is given as appendix A. 


1.- EQUATIONS OF SLENDER-WING THEORY FOR LIFTING WINGS 


In the slender-wing concept the spanwise extension of the wing is 
considered small compared with the length, and in the neighborhood of 
the wing the flow field assumes a two-dimensional character which for 
lifting wings is determined by solutions of the two-dimensional Laplace 
equation in planes transverse to the stream direction. The solutions 
for lifting wings are applicable throughout the Mach number range since 
the slenderness approximation causes the Mach number to vanish. Although 
the character of the solutions Is two-dimensional, the three-dimensional 
aspects of the wing problem enter through the expansion or contraction 
of the wing in the transverse plane and through the irrot at tonality con- 
ditions. The Kutta condition is required to render solutions of wing 
problems unique whenever the trailing edge is of the subsonic type, and 
in the slender-wing approximation this condition is imposed on all swept 
trailing edges. The character of the flow in a plane normal to the 
stream direction depends upon the admissible behavior of the velocity 
at the edges of the wing. Linearized theory requires that the velocity 
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velocity in the direction of the x-axis, the solution can be constructed 
from the complex velocity V(xj£) » v(xjy, z) - iw(xjy, z) where £ = y + iz 
is the complex variable and x occurs only as a parameter in the velocity. 
The quantities u> v, and w are the disturbance velocities in the x-, 
y-, and z-directions, respectively. 

In each of the three regions (sketches A and B) the solution can be 
reduced to a mixed-boundary- value problem considered by Cheng and Rott \ 
(ref. 7) } namely, the evaluation of a complex function in a half -space 
when the real and imaginary parts are prescribed alternately along seg- 
ments of a line. With the mean camber surface defined by z - Z(x,y) 
and with the subscript o used to denote values on z = +0, the real 
and imaginary parts of V(x}£) on z = +0 are alternately given in 

region I since w Q - — on the wing and v 0 is zero off the wing. In 

dx 

regions II and III, however, the real and imaginary parts of V(x;£) 
are not alternately prescribed on z = +0 since vo is not specified 
on the trailing vortex sheet. Differentiation of V(x;£) with respect 
to x gives a function with the desired properties since 

V x (x;£) = v x - iw x 

and by the irrotationality conditions, 


V x (x;5) = u y - iw x 

Thus the real and imaginary parts of V x (x;£) are prescribed alternately 
on z = +0 since vanishes everywhere off the wing, because the 

pressure coefficient must vanish there, and w^ is known on the wing. 

The prescribed values of the velocity components and their derivatives 
on 2=4-0 together with the requirement of vanishing velocity at infin- 
ity and the specification of the nature of the velocity at the edges are 
sufficient for the determination of the solution in each transverse plane. 


Solution of Mixed-Boundary-Value Problem 

The method of solution of the mixed- boundary- value problem for the 
upper half-space presented in reference J is based on the well-known 
integral relation 


V(£) = v - Iw 


= 1 
it 


s: 


vp(n) 

£ - T) 




d n (i.i) 
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for determining the complex function V(£) from a knowledge of the 
imaginary part on z = +0. Consider v 0 and v Q to be prescribed 
alternately along segments of z = +0. If a complex function r(£) 
can be determined, which is real on the segments where w Q is prescribed 
and imaginary on the segments where v 0 is prescribed, then the imaginary 
part of the function V(£) /r(^) is known everywhere on z = +0 and from 
equation (l.l). 


v<$) _ 1 

r i i 

[v(£)l 

r(5) * J 

-co i - n 

r(£)_ 


dn + p(s) 


z=+0 


( 1 . 2 ) 


where I denotes the imaginary part. 

The integral term in equation (1.2) does not provide a unique solu- 
tion since a function P(£) which is real on z = +0 can be added with- 
out altering the prescribed conditions. Cheng and Rott (ref. 7) have 
shown that, with w Q prescribed on the N segments ranging from b^ 

to b]£, the function r(£) is of the form 

4 “ /, xt^/, -\t +S * 

r ( ^ ) - i H - bjjJ - bj^J 


where m^. and are positive or negative integers. The order of the 

admissible singularities at the edges dictates the. values of m^ and 

for a particular problem. The function P(£) is a polynomial, in £, the 
degree of which is established by the requirement that the velocity vanish 
at infinity. The solution for the lower half-space is known by symmetry. 


Solution in Region I 

The application of equation (1.2) to region I gives the usual form 
of the slender-wing solution. From the boundary conditions on z = +0 
in region I, the real part of V(x;£) is zero for | y j > l and the 

imaginary part is prescribed for |y | < 1. Square-root singularities 
are permitted in the velocity at the edges, and the function r(£) is 
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. The requirement of vanishing velocity at infinity dictates 

Vs 2 - i s 

that P(xj£) is at most a function of x , and since there is no net 
circulation in the crossflow plane, this quantity is zero. From equa- 
tion (1.2), 


V(x*g) 




w 0 (x;ti) 


\JZH 

i - n 


drj 


(1-5) 


The velocity component v 0 (x;y) is given by equation (1.3) for 
-l ^ y S l on z = +0 and the velocity potential on the wing is obtained 
by integration as 


<Po(x;y) 


r y 

J v 0 (x,*n) 


dr) 


_ v* 2 y 2 


L 


v(x;ti) dt) 

1 \jl 2 - r) 2 (y - tj) 


(1.4) 


where v is given by the Indefinite integral v(x;y) = J' w 0 (x;ti) dTj 

since a constant can be added to v without altering the potential. 

The function w D (xjy) for y > l is required to determine the 
solution in region II for c < op and from equation (l. 3), 


i pl \Ji 2 _ -2 

v 0 (x;y) = = / v 0 (xjti)-J 1- di) (y > J) (3—5) 

. t a J-1 r-ji 

*Where no confusion can arise, as in equation ( 1 . 5)3 the argument 
of l (and also of the functions t and f introduced subsequently) is 
not explicitly indicated in the remainder of the report. The argument is 
either x or the streamwise variable of integration f- . The particular 
quantity in question is clear from the context and is stated explicitly 
for the final equations. 
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Solution in Regions II and III 


The solution in region II for c > c^ is structurally similar to 

that in region I and is easily evaluated. The determination of the pres- 
sure in region II for c < c^_ as well as the solution in region III is 

most readily accomplished by considering separately the cases of symmetric 
and antisymmetric spanwise loading. For these cases, only the function 
V x (x; 0 is evaluated in this section and the pressure is determined in 

sections 2 and 3* 


Solution in region II for c > C]_.- The trailing vortex sheet lies 

y[ > * 


in the region 

w °x 
for 


in region II for c > c^ ( sketch A) . The quantity 


Is known from the boundary data for -Z S y S l and Uoy is zero 
lyj > Z. The Kutta condition must be satisfied on y = ±Zj conse- 


quently, the velocity is regular and V x has a square-root singularity 


on ±1. The appropriate function r(£) 
tion (1.2), 


is — - , and from equa- 

\js 2 - * 2 


V x (x;£) 




£ - H 


( 1 . 6 ) 


The function P(x;£) is at most a function of x, and it can be 
shown readily that this function must be zero to satisfy the Kutta con- 
dition. The quantity' Uoy is given by the real part of equation (l.6) 

on z = 0, and since uo is zero off the wing,- the disturbance velocity 
on the wing is determined by integration from -Z to y as 


Uo(x;y) 



vq x (x;t)) 7] 2 


Interchanging the order of integration and integrating by parts gives 
the streamwise disturbance velocity as 
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uo(x;y) 




L 


v x (x;n) <hi 


i 


i 2 - n 2 (y - n) 


(i.T) 


where v x = J di). Comparison of equations (1.4) and (1.7) shows 

that uo(xjy) in region II for c > is given by the same integral 
expression as cp 0 (x;y) in region I with v(x;y) replaced by v x (x;y). 

In region II for c > c-j_ the slope of the edge y = l(x) Is less than 
zero and this region of the wing supports load provided w Qx does not 
vanish everywhere in region II. For the flat-plate wing, Wo x vanishes 

and sections downstream of the maximum span support no load. This is 
the result obtained by Jones (ref. l) . 

The function V x (x;£) in region II for c < c^.- In region II for 
c < c-^, the trailing vortex sheet lies between the two wing panels 
(sketch B) . The quantity Uo y is zero for |y[ > l and jy [ < t, and 
w 0x i3 known on the wing. Square-root singularities are admissible in 
the velocity on the leading edge; consequently, V x (x;£) has a three- 

halves -power singularity there while the velocity Is regular on the 
trailing edge and V x (x;£) has a square-root singularity. The appro- 
priate function r(£) is - , and from equation (1.2), 

V x (x; 5) in region II for c < c-j_ is 


V x (x;£) 




f i VoxUn) (j a -Tl 2 )^T] 2 -t 2 wc^Cxjii) (l 2 -n 2 )^ 2 -t 2 ^ + P(xj£) 


Cl. 8) 


where the coefficients in the polynomial are functions of x. 

The function V x (x;£) in region III.- The solution in region III 
differs from that' in region II for c < only in that the velocity is 
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regular on all edges in the former since the Kutta condition must be 
satisfied on both l and t. The function r(£) then is 

3- and from equation (l.2), 

\J(? 2 - i 2 )(t 2 - 1 2 ) 

«, D . t ■ g; fc * H <i - w 

2.- WINGS WITH SYMMETRIC S PANWISE LOADING 

The expressions for the potential in region I and the Jtaturbance 
velocity uo In region II for e > Cl can be specialized for either 

symmetric or antisymmetric spanvise loading end the 

pressures are obtained directly from these quantities. In region II 
for c < cn and in region III it proves convenient to make use of t 

spanvise symmetry conditions at the outset to J ^^J^Sped 

sures. The surface pressures for symmetric spanvise loading are devexopea 

in this section and those for antisymmetric loading are presented 
section 3* 

For symmetric spanvise loading v Q (xjy) = wo(x,-y), and the surf 
pressure coefficient in region I is determined from equation (1-4) as 


(2.1) 


r l 1+ iEf=2T z 

r Pm 2 / v ( x ' T| ) log - d^ + f WoCxjti)-^: 

cp = -^o x - * w °x (x ’ n)log ^ p-r ^ 

U 2 -n 2 


vhere the prime on l denotes the derivative with respect to x. In 
region II for c > from equation (.l-7;> — 

Op - -2=o-| [ v 0x (x ;1 )lo g 1 " M =£ «. ^ 

Jo 1 - Y, 2 _ n 2 


4 
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Prom these equations the pressure is seen to be discontinuous at the 
boundary of regions I and II for c > c 1 if l'(c 1 ) is different from 

zero. 


Solution in Region II for c < 

Pressure.- For symmetric spanwise loading the integration of equa- 
tion (1.6) can be carried out over the right wing panel to give 


V x (x;0 = 




p(xjg)- 2^ v 0 x (x}n).(i 2 -n 2 ) 2 jn^t^ 


(2. 


The symmetry condition requires that the polynomial be an even function 
of £, and for the velocity to vanish at infinity, P(xj£) = A + 
where A and B are functions of x. The quantity Uoy is given by 

equation (2.3) for t ^ y ^ l, and since Uq is zero off the wing, the 
disturbance velocity on the wing for t = y ^ l is 

py 

uo(x;y) = j- ^ dT l 


and after some reduction 


uo(xjy) = 


* l 


{i 2 - 1 2 ) 


( A+BZ2)e( 0, k) - (A+Bt2)F( f, k)+( A+Bl2)£fefl£ 


l h 2 -y 2 


-G(x;y) 

(2.4) 


where F(0,k) and E(0,k) are incomplete elliptic integrals of the 


I 2 - y 2 

777 


first and second kind, respectively, sin 


the modulus k 
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is given by k? = 1 and 


G(x;y) - f £ 




y (i 2 - n 2 ) 3/2 f 2 - 12 


v 0 x (x;n) (* 2 - tj 2 ) ^tj 2 - t 2 — — 


ti 2 - TJ 2 


(2.5) 


The sign denotes the finite part of an integral. Satisfying the Kutta 
condition at y = t gives the streamwise disturbance velocity as 


uo(xjy) 4 - + - - BZ 


itz(z 2 - t 2 ) 


Z(tf,k) + 2 , 


- t c 


-y 2 


+ <K*t) 


G(x;y) 


( 2 . 6 ) 


where Z(0,k) = E($,k) - ^ F(0,k) is the Jacobian zeta function. 

K 

The functions G(x;y) and G(x;t) are evaluated in appendix B. 
These functions can be reduced to a single integral involving elliptic 
integrals and related functions. The combination of these two functions 
which occurs in equation (2.6) can be expressed in terms of the quantity 
in brackets in equation (2.6) and a theta function. With the use of 
equation (B4) the pressure coefficient in region II for c < c^ then 

becomes 


C p - 


— 

-2S(x) l ' Z(0,k) + + * f t dT l 


(2.7) 


where the unknown quantity S(x) which is related to A and B by 
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S(x) l ' - 


nl{l 2 - t 2 ) _ 


2 / ^oxU^) ^(l 2 - 1} 2 )(t] 2 - t 2 j T] dTi - (A + Bl2) 


is introduced to facilitate the analysis. The function fi(xjy,Ti) is 
expressed in terms of an infinite series by equation (B3) . Equation (2-7) 
i3 equivalent to that obtained for the pressure coefficient for the flat- 
plate wing in reference 5 with the addition of the integral involving 
ft(xjy,Ti) and the integral in the definition of S(x) . 

Integral equation for S(x) .~ The Kutta condition gives one relation 
between A and B and an additional relation is required to evaluate 
these functions, or equivalently, the function S(x) . Since the solution 
has been constructed from Wo x whereas the surface boundary condition is 

in terms of w 0 , the second relation is that the solution must give the 
prescribed v Q on the wing. The function S(x) is determined by a 
procedure essentially the same as that used by Mirels (ref. 5). First, 
w Q on the wing is expressed as 


w Q (x;y) * J w O|(£?y) U 


where it is essential that y be greater than Z(c) for the path of 
integration, since for y < Z( c), the relation is an identity and gives 
no condition on S(x) . Separating the integration into the three parts 
-oo to c, c to Z’l(y), and Z”^(y) to x gives the equation 


WolV^yby) - Wo(cjy) 


l 


^ _ 1 (y) 


w 0 e (l;y) 


( 2 . 8 ) 


where w 0 ^Z -1 (y);yj is the prescribed value of w 0 along the leading 
edge and w Q (c;y) is the value of the downwash off the wing in region I 
and is given by equation (1.5) with x = c. 

Equation (2.8) provides an integral equation for the determination 
of S(x) but involves both the unknown function S(x) and dS/dx; the 
function ds/dx arises in the course of evaluating the finite part of 



NACA TIT k2k2 


13 


the integral. An equation in terms of S(x) only is obtained by inte- 
grating equation (2.8) from y = l(c) to y= l(x) . Then, 





w 0 ^(iiy) fly 


(2-9) 


where the order of integration of the double integral has been inter- 
changed. The region of 
integration of the right- 
hand member of equa- 
tion (2.9) is shown 
shaded in sketch C. 


The expression for 
appearing in the 


W 0 £ 

integral of equation (2.9) 
is obtained from equa- 
tion (2.3) for y > l. 
Satisfying the Kutta con- 
dition from equation (2.k) 
gives a relation between 
the functions A and B, 
and the polynomial can be 
expressed in the form 



P(x}£) s A + B£ 2 = -nl(t 2 - l 2 ) 


A + Bl 2 


nl(l 2 - 



* 2 


G(x;t) 

K 


where K and E are the complete elliptic integrals of the first and 
second kind, respectively, with modulus k. Then, when the integration 
is performed, the inner integral in equation (2.9) becomes 



<2y = s(ljKx)) 


+ > G(Ut) + ^ 


Ao(e,k) 


X(x) 


>/l2(x) - l*\fl2(x) - t2 


( 2 . 10 ) 
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where Ao(0^k) is Henman's lambda function (ref. 10), k 1 is the 
complementary modulus, and 


k' = \jl - k 2 = | 


sin 8 = 


l 2 (x) - l 2 
Z 2 (x) - t 2 


The function G(xjt) is given by equation (2-5) with y = t and 


s(ftJ*(x)) 


( y 2 _ l 2 ^/ 2 ^ 2^2 1 ' y 2- n 2 


The functions G(|jt) and g(|}2(x)) are evaluated in appendix B. 
These functions can be reduced to a single integral with the integrand 
expressed in terms of elliptic integrals and related functions. The 
combination of the two terms which arises in equation (2.10) can be 
expressed in terms of the quantity in brackets in equation (2.10), a 
theta function, and elliptic integrals. Prom equations (2.10) and (B6), 

f I(x) *o.(tiy) dy = H(S;l(x)) - S(t)I'(|)[~f — - SjL*i 

Jl ‘ |_ 2 _ K fl 2 (x)-I 2 ^(x)-t 2 


where 


F(e,k T )F(oo,k) r(SiTj) 


In equation (2.13), K' Is the complete elliptic integral of the first 
kind with the complementary modulus k' and sin ai = 

9 
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function r(|;Tj) is evaluated from, the series given by equation (B9) • 

From equations (2.9) and (2.12) the final form of the integral 
equation for S(x) is 



it -Ao(®>^) t k 2 Z l(x) 

j u J 

* c 

2 K \/z 2 (x) - l 2 ^ 2 (x) - t2 


d| (2.11) 


where the left-hand member is 


M(x) = J' w Q (c;y) - v 0 (l~Hy) »yj dy + J' H(|;Z(x)) d| ( 2 . 15 ) 


1(c) 


w 0 (c;y) is given by equation ( 1 . 5 ) with x = c, and 

is the prescribed value of w Q along the leading edge. 

In equations (2.13) and (2.ll) l, t, and k are functions of the 
variable of integration | where the argument is not explicitly indi- 
cated. The kernel of equation (2.14) is equivalent to that given by 
Mirels for the flat-plate wing at angle of attack a. For the flat 
plate, w^ vanishes; hence, the function H(|; l(x) ) and the integral 

in equation ( 2 . 5 ) are zero and M(x) = a jj l 2 (x) - Z 2 (c) . 

Equations ( 2 . 7 ) and (2.14) together determine the pressure in 
region II for c < c-j_. Equation ( 2 . 14) can be inverted by numerical 

methods to obtain S(| ) by a method similar to that of reference 5* 

(See appendix C.) The Jacobian zeta function and the Heuman lambda 
function are tabulated functions (refs. 8 , 9, and 10, for example). The 
functions &(x;y,T}) and r(|;T}) which arise in the equations for the 
pressure coefficient and M(x), respectively, are determined from an 
infinite series. These series converge rapidly, however, and the evalua- 
tion of fi(x;y,T)) and r(§;T]) offers no serious difficulty. 

Calculation of the pressure for x near c . - The integral equa- 
tion (2.14) and the expression for the pressure coefficient (eq. ( 2 . 7 )) 
can be approximated to give a somewhat simplified as well as useful 
solution for x near c. In this neighborhood, the kernel of the inte- 
gral equation can be approximated and the inversion for S(x) can be 
obtained in a closed form. The integral equation is approximated by 


The quantity 

w 0 (^ _ 1 (y);y) 
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which is equivalent to that given in reference 5 . The inner integral is 
the logarithmic integral (tabulated in ref. 11 , for example) and 
cL*t 

7 = — • The quantity in brackets is independent of x, and the 

cU x=c 

inversion of this integral equation of the Abel type is 


S(x)Z(x)Z '(x) = i 


ill 

iJSSclKxJl-OO f 4 

k [_ 7 «0 log n 31 dx J c 


M(|)U' cL£ 


? (x) - Z‘ 


(2.16) 


where k and k' are functions of x. Frcxn equation (Bll) the func- 
tion H(£;Z(x)) for x near c is 


H(i}Z(x)) 





n Jn 2 - t 2 

lh 2 - T , 2 


+ Z(co,k) 


dT] 


(2.17) 


The function &(x;y,Tj) which arises in the evaluation of the pressure 
is approximated by the first term of the series representation as 


fl(x}Tj) 


a(x; tj) 


— log 
2 


— log 
2 


sin(a + t) 
sin(cr - t) 



(x near c) 


> 


(x = c) 


( 2 . 18 ) 


j 
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where 


^ _ jC_ k) 

2 K 


T= it F(gU) 

2 K 


It can be shown that S(c) depends on the c amb er distribution at 
x = c through the relation 


S(c) 



•w q (c;ti) cLu 
\/l 2 (c) - ti 2 


and since 

and (2.7) 
regions I 


Z(0,1) 



comparison of equations (2.1) 


shows that the pressure is continuous across the boundary of 
and II for c < c-j_ if X '(c) is continuous. 


Solution in Region III 

For symmetric spanwise loading equation (l-9) can be written as 


v x (*;0 = 




Hx; i) + 2 J w 0x (x;n) \Z( l2 -i a ) (n 2 -t 2 ) 


(2.19) 


The condition of symmetry and the requirement of vanishing velocity 
at infinity determine P(x;£) = A where A is a function of x only. 
The quantity Uq v is given by the real part of equation (2. 19 ) , and 

integration gives 


uo(x;y) 



^(i 2 -u a )(u 2 -t 2 ) 



vo x (x;ri)^(l 2 -n 2 ) (n 2 -t 2 )-l^| 



du 
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Satisfying the Kutta condition at y = t determines A. Interc han ging 
the order of integration and integrating gives the pressure coefficient 
in region III 


C P 


_4_ 

t a 



w 0x (x;n) 



F(9,k) 

K 



n 


dT] 


4 

it 



w o x (x}Ti)n(x;y,T]) dT) 


( 2 . 20 ) 


where Il(e,a 2 ,k) is the elliptic integral of the third kind with the 
l 2 - t 2 

parameter a 2 = — — -j the last form of equation (2.20) follows from 

1 2 „2 

- n 

the substitution of equations (B2) for the elliptic integrals. From 
equation (2.20) the pressure coefficient in region III is zero only if 
w 0x is zero. Comparison of equation (2.20) with equations (2.1) 

and (2.7) shows that the pressure across the boundary of regions II 
and III is discontinuous for c < if l'(c 1 ) / 0 and is continuous 

for c > c-^. Thus, for either c > c^ or c < c-]_, the surface pressure 

Is discontinuous across the plane normal to the Btream direction which 
passes through the wing maximum span unless l * Is continuous and zero 
at x = c^. Further, the surface pressure is continuous across the plane 

passing through the root chord if 2'(c) is continuous. 


Lift 

The lift can be determined from the asymptotic form of V(xj£). 

For a wing, the function V(xj £) for large values of £ is of the form 


V(x;£) = - ^mamS"^ 1 ) (2.21) 

m=l 

From reference 5 or 12 the lift carried by the wing upstream of x is 
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L = I^itq^aq) 

where Is the stream dynamic pressure. 

The coefficient a-j_ in region I is determined from equation (l.J)} 
thus L becomes 


L - 'Ho J V 0 (x}T])\Jl 2 - T\'‘ 


dq 


( 2 . 22 ) 


In region II for c > cq 
equation (1.6), and dL/dx is 


the function daq/dx 
then determined as 


is obtained from 


dL 

dx 




di) 


(2.23) 


Similarly, the function dL/dx in region II for c < c 1 is determined 
from equation (2.3) as 


dL , 

— = ^t<lj 


S(x)l 




w 0x (xjTi)Z(a),k) dq 


( 2 . 2k) 


From equation (2.19), dL/dx in region III is 


dL 

dx 



w 0x (x;Ti)Z(cD,k.) dr) 


(2.25) 


and the total lift is obtained by integration. 

The lift supported by the part of the wing in region I is dependent 
only upon Wo(cjy) (or w 0 (cq;y) for c > cq) . In regions II and HI, 

however, the contribution to the lift depends upon the camber distribution 
throughout these regions. 
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Wings With Twist 

The slender-wing solutions assume their simplest form, for the flat 
plate and. for wings with twist but no camber. For these conditions the 
wing supports load only where Z' > 0. The flat -plate solution has been 
given in references 5 and 6 , and the solution for twist proportional to 
y2 is developed here to illustrate the use of the equations. 

Let w 0 (x;y) = — y2 where a is a constant. Then from equa- 

‘(c) 

tion (2.1), the surface pressure coefficient for x = c is 


C p = 


-29o x - 


Z 2 (c) 


Inl- 




and from equation (2.22) the lift carried by the part of the wing upstream 
of x in region I is 


L « 



The loading is zero in regions II and III except for region II with 
c < cp. The solution in region II with c < cp requires the evaluation 

of the function S(x). Only the terms w 0 (c;y) and w 0 (z _ ^(y) ;y) con- 
tribute to the left-hand member of the integral equation for S(x) ; from 
equation (1.5), 


w 0 (c;y) 



2y3 - yZ 2 (c) - 2y 2 


^y 2 - Z 2 (c) 


and w 0 (z -1 (y) jy) = — y 2 . Then from equation (2.15)., the left-hand 

Z 2 (c) 

member of the integral equation is 

M(x) = - — [ 2 Z 2 + Z2(c)l \/z 2 - Z 2 (c) 

Sl 2 (c) L ' 
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The inversion of equation (2.l4) for x near 
tion (2.16) as 


c 


is given by equa- 


S(x) 




kj_ 

a 

4 

f 4 du ... * 2 (x) 

2k 

7 v. 

/ 0 l0 S V- l 2( c ) 


( 2 . 


and S(c) = - — . 


2 

for c < x < c^ 


From equation (2.7) the surface pressure coefficient 
is 



-2S(x) 1 1 


z(?U) 


+ 


£ 

1 



and at x = c. 


Cp “ 


- 7 


From equation (2.25) the lift per unit length for c < x < cj_ is 


and at x = c, 

g= 2x^11' 

Comparison with the expressions for the lift and the pressure coef- 
ficient in region I shows that dl/dx and Cp are continuous across 

the boundary of regions I and II if l 1 is continuous at x = c . 
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The function S(x)/a computed from equation (2.26) is compared in 
figure 1 with that calculated by the numerical method given in appendix C 
for a wing with straight edges (l = m^x and t = mfc(x-c)) for 7=1, 

2, and 3* The interval for the numerical method was taken as 0.05, and 
Simpson's rule was used to evaluate the integrals. The results show 
that the approximate inversion provides a very satisfactory solution 
even when k differs appreciably from unity. 


3.- WINGS WITH ANTISYMMETRIC SPANWISE LOADING 


The pressures in region I and region II for c > c^ for antisym- 
metric spanwlse loading follow directly from equations (1.4) and (1.7). 
The pressure in region II for c < c^_ and the pressure in region III 

can be determined in a manner similar to that for the case of symmetric 
loading. 

For antisymmetric spanwise loading, w 0 (x}y) = -w D (xj-y) and the 

surface pressure coefficient in region I is determined from equation (1.4) 
as 


C p = 


-aPox = | / ^(xjtiJlog 


1 

>. 

1 

cva 

rJ 



,2 2 

1 - l] 

l -y 

y 

,2 2 
l -ip 


dT] + 


^ f ' 


VoUinh dq 




(5-1) 


In region II for c > c^_ from equation (1.7), 


Cp - 


- 2u ° = i 


Wo^XiT^lOg 


l + ±' 11 


- y 2 


l-ll 


y yz 2 - t ] 2 

1^7 


dq 


(5-2) 


y w 1 2 


Thus the pressure is discontinuous at x = cj_, the boundary of regions I 
and II for c > c^, if 2'( cq) is different from zero. 
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Solution in Region II for c < cq_ 

Pressure in region II .- For antisymmetric span-wise loading, equa- 
tion (1.6) becomes 


v x (xjt;) = 


-i 


«(t 2 -i 2 )\ |^ L 


^ ~ 2 + 


( 3 - 3 ) 


The symmetry condition requires that the polynomial be an odd function 
of £, and since the velocity must vanish at infinity, P(x$ 5) = + B£3 

where the coefficients A and B are functions of x only. The quan- 
tity UQy is given by the real part of equation (3- 3) , and since Uq 

is zero off the wing, the disturbance velocity on the right wing panel 
is 


uo(xjy) = jf Uo^(x;p) dp. 


jy 2 -. t2 + = ! B ±n- 1 J l2 - y2 . + -i- 2 ^" t2 
l 2 - y 2 * \ fl 2 -t 2 l 2 -t~[l“-y“ 


-•f f Vn ( x ; ti ) ( i 2 -^ 2 ) 2 jjn 2 -t 


2 

V- 2 -t\ 2 


With an interc han ge of the order of integration and after integration 
on p, it can be seen that the Kutta condition requires that B equal 
zero. The linearized expression for the pressure coefficient on the 
right wing panel then becomes 
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Cp = 


-2uo = 


2 
j ( 


a 

■w 0 x (x;T])log 




(y2-t 2 ) 

¥7) 

f_2 +2'| 

(n -t , 

j 

N 

(y 2 -* 2 ) 

If 

(i 2 -* 2 ) 

U 2 -t 2 ) 


2Q(x)U 



( 5 - 4 ) 


and is of the opposite sign on the left panel. The function Q(x), which 
is related to the coefficient A by 


Q (x)U' 



\j{ l2 - n 2 ) (n 2 - t2 ) dT i 


is introduced for convenience. 

Integral equation for Q(x) . - The function Q(x) can be determined 
from an integral relation in a manner similar to that -used to evaluate 
S(x) . The procedure is to integrate w 0x over the shaded region in 

sketch C. Then the integral relation is 


.l(x) 


IZ K*- ,;L (y);y) - w 0 (c;y) 


fly 


pX pl(x) 

= J d ljf^ v 0 g (fc;y) dy 


( 5 - 5 ) 

where w 0 ^Z -d (y);y) is the prescribed value of w Q on the leading edge 
and w Q (c;y) is given by equation ( 1 . 5 ) with x = c. Substituting w D ^ 

from equation (3*5) into equation (3- 5) and integrating on y gives the 
integral equation for Q(x) 


N(x) 





'; 2 (x) - t 2 
l 2 (x) - z 2 




( 5 . 6 ) 
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where the left-hand member is 


N(x) = J v c (c;y) - v 0 ^ _1 (y);y) dy + J h(|jl(x)) d| (3-7) 


1(c) L 


The quantity w 0 (cjy) is given by equation ( 1 . 5 ) "with x = c, 
w 0 ^Z"^(y);y) is the prescribed value of w Q Along the leading edge, 
and 


/ [l 2 (x) - 12] 

(n 2 - t 2 ) 

/ _Z 2 (x) - t 2 ] 

(i 2 - , 2 ) 


dt] (5-8) 


In equations (3*^)j (3-6), and (3-8) the argument of l and t is | 
except where indicated explicitly. 


Mirels ' solution of the flat-plate rolling wing is a particular 
example of these equations. For the rolling wing w D = -tr^y, where o^. 

is the angular velocity, and since w Q is independent of x, the func- 


tion h(|jl(x)) 


vanishes and N(x) 


be comes 


| a^l^l 2 - l 2 (c). 


The numerical solution of equation (3-8) is readily determined by 
the method of reference 3 once the function N(x) is evaluated. (See 
appendix C.) For x near c, t is small compared with l and can 
be neglected, and the solution is given in closed form as 


Q(x)Z(x)Z*(x) 


2 _d_ 

3 1 dx 



nQ)z’ as 
V* 2 (x) - i 2 


(5.9) 


The function Q(c) can be shown to depend upon the camber distribution 
at x = c through the relation 


Q(c) 




w 0 (c;t])ti dt] 
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Solution in Region III 

For antisymmetric spanwise loading, equation (1.9) can be -written 


as 


v x (x;5) = 


W(s 2 - i2 )(s 2 - t2 )L 


^o x (xir,)\J(l 2 -n 2 )(Tl 2 -t £ ) * A? 


(5-10) 


where the polynomial is at most A£ from the symmetry conditions and 
the requirement of vanishing velocity at infinity. For t £ y ^ Z on 

z = +0, equation (3*10) gives Uq , and Uq on the right wing panel is 

i7 

obtained by integration from 2 to y as 


uo(x;y) = - | f 

J V 


H dp 

y \/(i 2 -, 2 )(, 2 -t 2 ) Jt 


f t V*” 1 * \|^-n 2 )(l 2 -t 2 )^ 


where A is determined as zero by the Kutta condition at y = t. With 
the use of the linearized pressure relation and after an interchange of 
the order of integration the pressure coefficient for t S* y ^ l becomes 



(3.H) 


and is of opposite sign on the left panel. Comparison of equation (3- 11) 
with equations (3-1) and (3*4) shows that the pressure across the boundary 
of regions II and III is discontinuous for c < c-^ if 2'(c]_) ^ 0 and 

is continuous for c > c-]_. 
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The distinguishing feature of slender wings with cambered surfaces 
as compared with flat-plate wings, for either spanwise symmetric or anti- 
symmetric loading, is that cambered surfaces carry load downstream of 
the position of maximum span. The pressure is discontinuous across the 
plane which passes through the wing maximum span unless l ' is continuous 
and zero at x = c^_. Further, the pressure is continuous across the plane 

through the root chord unless l 1 (c) is discontinuous. 


Rolling Moment 

The moment about the x-axis due to antisymmetric spanwise loading 
can be evaluated in a manner similar to that for the lift (ref. 5)- With 
the asymptotic form of the complex velocity given by 


V(x;0 = - ^ 
m=l 


the moment 


about the x-axis is 




The contribution to the moment in region I for the part of the wing 
upstream of station x, with the use of equation (l.J), is 


My 


= v o ( xjt ])\j 


r\ dT) 


In region II for c > c^, the quantity ds^/dx. 
asymptotic form of equation (l.6)j then, 


( 3 . 12 ) 


is determined from the 
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and in region II for c < cq, da^dx is evaluated from equation (3.3) 
and 


£r = atq*. (l 2 -t 2 ) Q(x)U> - 2_jT Vo x (xiTi) \J(l 2 -r} 2 ) (r} 2 -t 2 ) d^ 


(3-14) 


Similarly, in region III 


dx 


= -Ho J v o x (x;t])^1 2 - T \^( t \ 2 - t 2 ) dTl 


(3.15) 


The contribution to the rolling moment for the whole of region I depends 
only upon w G (c;y) (or w Q (c 1 ;y) for c > c x ), whereas the contribu- 
tion to the moment in regions II and III depends upon the distribution 
of camber throughout these regions. 


4.- CONTROL- SURFACE SOLUTIONS 


The evaluation of the pressure for cambered surfaces involves con- 
siderable calculation since, in general, the various integrals must be 
evaluated by numerical methods. The additional load due to the deflec- 
tion of a control surface, however. Is readily determined from the gen- 
eral equations since many of the integrals can be evaluated analytically. 
Aside from the overall slenderness conditions required for the wing, the 
control surfaces should also satisfy the same conditions. Solutions are 
given for both symmetric and antisymmetric spanwise loading and are 
denoted as flaps and ailerons, respectively. These expressions apply 
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to trailing-e&ge control surfaces with side edges such as those shown in 
sketches D and E; similar 
expressions for tip con- 
trols can be obtained by 
the same method. 



The leading edge of 
the control is defined by 
y = f(x) = m f ^x - x f ). 

On the right wing panel, 
f' > 0 for the controls 
of sketch D and f 1 < 0 
for the controls of 
sketch E. The coordi- 
nates of the upstream and 
downstream leading-edge 
corners are ^X]_,y^ and 

( x 2>y2jj respectively, 

\ 

and the control-surface 
deflection is 5 C and is 

positive for a downward 
deflection. 

The additional load 
due to the deflection of a 
control surface is obtained 
from the general equations 
for cambered surfaces by 
setting w Q = 0 on the 

wing and w Q = -S c on the 

control surface. The quantity w^ is zero everywhere except on the 



edge y = f(x) 
which contain 


where w Q is discontinuous. In each of the integrals 


w, 


^ as a factor, the quantity w Qx is interpreted as 


the component of a Dirac delta function, 
the form 


and all such integrals are of 


/ 


w o x (xJ T l)^( x J T l) 


dii = -6 c |f’|F(xjf) 


where the integration across the boundary y = f is in the positive 
y-directionj the value of the integral is zero if the integration does 
not extend across the edge y = f . The expressions for the loading apply 
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to the configurations of both sketches D and E, with the (+) sign 
applicable for f ' > 0 (sketch D) and the (-) sign applicable for 
f' < 0 (sketch E) . The equations for the loading are directly appli- 
cable for xi in region I and x 2 in region III. The additional loa ding 

is zero upstream of xp for any location of x^> and the solutions for 
x > x 2 , regardless of the location of x 2 , are obtained from the equa- 
tions presented by replacing f by y 2 and setting f' = 0. 


Flaps 

Region I .- The pressure is zero upstream of x-j_ since w Q = 0 on 
the wing, and from equation (2.1), 


C P 




From equation (2.22), the lifting force in region I is 


L - 


±^<2oo 8 C 


f 




+ l 2 



f 

l 



and the total lift in region I is obtained by setting x = c where 
c < c-j_ and x = c^ where c > c^. 

Region II, c > c-j_.- From equation (2.2), the pressure coefficient 
in region II for c > c-]_ Is 


Cp - - 



log 


1 + 


H 2 - 


1 - 
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The lift per unit length is evaluated from equation (2.23) as 




w 


l 2 - f 2 


and the total lift in this region is found by integration from to c. 

Region II, c < Cq_.- The pressure coefficient for c S x ^ c^ is 
given by equation (2.7) as 


Cp - 


2S(x) l ' 

ztf,*) h- y|/»® - *? 

_ 

f 1 


1 1 l 2 - y 2 . 

it 



n(xjy,f) 


where the function fi(x;y,f) is defined in appendix B with tj replaced 
by f. 

The function S(x) is evaluated from equation (2. 14). The quantity 
w Q (z - l(y) ;y) is zero and the left-hand member of the integral equa- 
tion ( 2 . 15 ) is evaluated as 


M W . ± ^ L> nt, - a ^ ■ ,2<<o ■ “"' x ■ 




f(e) . / - X g (c) 

1 Vi 2 (=) - 


tan -1 


yx. i 2 - iH c) 


1 u,2, 


i^Cc) - y x c 


- J hUjI(x)) 




(*• 


where H(|;Z(x)) is determined from equation (2.13) as 


H(|;Z(x)) 


= -5c 


f 



F(9,k l )F((ju,k) 

K'K 


r(S;f) 

K 1 


and the function r(|;f) is given in appendix B with tj replaced by f. 
If X£ lies in region II, the upper limit in the integral in equa- 
tion (4.1) is X2 for x > X 2 - 
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The leading terms in the function M(x) are 

whereas the contribution from the integral term is 
and for x near c. 


o(l 2 - Z 2 (c)) 1 / 2 
o(l 2 - l 2 (c)) 3/2, 


r 


M(x) » ± 


26 f 


jt 


4 _i f(c) 4 _i y l 

BiXL 1 TT^T ~ sin 1 T(c) 


\ji 2 - i 2 (c) + o(i 2 - iHc)f 2 


The term o(l 2 - l 2 (c))^/ contributes to dS/dx but not to S(x) at 
x = c. Higher order terms can be retained in the approximation to M(x) 
to improve the representation of S(x). An alternate procedure which 
facilitates the calculations is to compute the function M(x) from 
equation (4.1) and then represent it in the form 


M(x) 



sin 


-1 f(c) 

1(c) 








(4.2) 


With this representation the f -unction S(x) is evaluated from equa- 
tion (2.16) in the neighborhood of x = c as 


S(x) = Si 


sin” 1 - f ( c ) - sin -1 -~r 
1(c) 1(c) 


1 + 


-f 

7 Jo 


4 ,.j3p ) + 

log tx y 


m 

V !*_ 
L ? 2n 


rt=l 


(2n + l)i 

I 2 - l 2 (c) 

(n!) 2 

1H0) J 


(4-3) 


and at x = c, 


S(c) = ± 


26 r 


j, _i f(c) 4 _1 ^1 
sin x ■ — ■ ■ - sin x 


1(c) 


1(c) 



z 
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The function S(x) can be evaluated by the method of appendix C in 
regions where this approximation is not suitable. 

The lift per unit length is determined from equation (2.2^) as 

f'|z(co,k) 


region II for c < C]_ 


Region III .- The pressure coefficient in region III for x < Xg 
is evaluated from equation (2.20) as 

Cp = - ^^£|f ' |fi(x;y,f) 

and is zero for x > xp* The lift per unit length is obtained from 
equation (2.25) as 


dx 


S(x) Z 


•M)*? 


where sin co = 


IZ 2 - f 2 


t 2 - t 2 


, and the total lift i 


is obtained by integration from c to c-j_. 


g= 8q e0 5 c l|f'|z(a),k) 

with Tj replaced by f in the definition of to, and the total lift 
over region III is obtained by integration from C]_ or c, whichever 

is larger, to xg- 

From the foregoing equations and the modifications for x > X2> it 

can be seen that the additional loading due to the deflection of a flap 
is different from zero only for x^ < x < X£ or c^_, whichever is larger, 

regardless of the location of x^_ and xg- Further, the pressure and 
lift distributions are discontinuous at xj_ and Xg and at cj for 
x 2 > c i if Z 1 ( c x ) 7^ 0 . 
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Ailerons 


The equations for the additional loading due to the deflection of 
an aileron are developed in the same manner as those for the flap except 
that the basic equations are those for ant isyxmne trie spanwise loading. 

The solutions for the aileron are somewhat simpler, however, since they 
involve only elementary functions. The equations presented apply directly 
to the configurations of sketches D and E with X]_ in region I and xg 

in region III. The additional loading is zero upstream of x^_ for any 
location of xq, and the solutions for x > X£ for any location of xq 
are obtained from the equations presented by replacing f by yp and 
setting f' = 0. 

Region I .- The pressure coefficient is zero upstream of x^, and 
for x > xq, the pressure coefficient on the right wing panel is evalu- 
ated from equation (3-1) as 



From equation (3-12), the rolling 


moment is 


^x = 





and the total moment from region I is obtained by setting the argument 
of l and f equal to the smaller of c or c^. 

Region II, c > c^. - The pressure coefficient in region II for 
c > c-j_ on the right panel is evaluated from equation (3-4) as 
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2. + -\ 

fiLud 

_1£] 

12 . f 2 


,2 2 

i - 

l - y 

y il 

2 2 
l * - f* 


From, equation (5*13) > the rolling moment per unit length is 


“ Ws c |f’|f\(j 2 - f 2 


and the total contribution to the moment is obtained by integration on 
x from c-j_ to c. 


for 


Begion II, c < c-j_.- From equation the pressure coefficient 



where Q(x) is determined from the integral equation (5*6). The quan- 
tity v 0 (l~Hy) is zero, and from equation (3-7 ) j the left-hand 
member of the integral equation is 
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where 



h(|;2(x)) is evaluated from equation (5.8) as 


cL| 

(4.4) 


. 26 c I I , - 2^](f^ - t^) 

, ,4, . 


If X£ lies in region II for c < c-p the upper limit in equation (4.4) 
is X2 for x > X2* The dominant terms in equation (4.4) are 

/ o ,\l/2 


I 2 2, ,W 2 

012 - l {c) J for x near c and the contribution from the 

integral is 0^2^ - iHc)Y /Z ; only the terms of 0 (l 2 - l 2 (c)) l/2 
contribute to Q( c) but the higher order terms contribute to dQ/dx 


at x = c. 


The function Q(x) for x near c is conveniently evaluated by 
computing N(x) from equation (4.4) and then representing it in the 
form 
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The quantity in braces is the term 0 (l2 - 12(c)) 1/2 obtained by 
expanding equation (4.4). From equation (3*9) then. 


Q(x) l = 1± 

31 



f 2 (c . )„ 


yi 


|/l 2 (c)-f 2 (c) ^(cj-y^ 


- 2 . 


m 

I 

n?=l 


an 


( 2n+l ) l 1 


2 

2 2 -1 2 (c) 

1(c) 

_l 2 (c) 


]jl 2 ( c)-f 2 (c) \jl 2 (c)-y^ 




(4.6) 


and 


Q(c)Z(c) 



\jl 2 (c.) - f 2 (c) 


The function Q(x) can be evaluated by the method of appendix C in 
regions where equation (4.6) does not provide a suitable approximation. 

The moment per unit length is evaluated from equation (3*14) as 


dx 





(^• 7 ) 


and integration on x from c to c^ gives the contribution to the 
moment from the part of the wing in region H. 

Region III .- From equation (3-ll)> the pressure coefficient in 
region III for t ^ y £ 1 is 
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and is zero for x > x 2 . Evaluating the integral in equation ( 5 - 15 ) 
gives the rolling moment per unit length as 

j—- Ho&c|f f | \]{} 2 - f 2 )(f 2 - t 2 ) 


and the moment from the part of the wing in region III is found by inte- 
gration on x from c or c^, whichever is larger, to x 2 . 

As in the case of the loading due to a flap, the additional loading 
due to the deflection of an aileron is different from zero only for 
x x < x < x 2 or c-l, whichever is larger, regardless of the location of 

Similarly, the pressure and lift distributions are discontinuous 
at x-l and x 2 and at e 1 for x 2 > c-|_ if 4 °* 


Example . - Calculations have been made for the swept -wing configura- 
tion with an aileron 
shown in sketch F. The 


leading and trailing 
edges of the right wing 
panel are defined by 
2 = m^x and 

t = m^(x - c), respec- 
tively, and the relevant 
constants are 

m^ = 0.5 mt = 1.5 

Elf =* 1 c l/ c = 1*2 

x l/ c = 0.9 x 2 /c = 1.25 
Xf/c = O.75 s/c = 0.6 
y-j/c = 0.15 y 2 /c = 0.5 
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QH • 

The quantity — calculated by the numerical method of appendix C 

Sc 

is compared in figure 2 with that calculated from equation (4.6) with 
m = 3. An Interval of 0.02 was used in the numerical method. The coef- 
ficients were evaluated so that equation (4.5) gave the correct value of 
N(x) at x = 1.02, 1.06, and 1.20. The results show that the approxi- 
mate equation for Q(x) provides a very satisfactory solution even when 
t/l is not small compared to unity. For example, at x/c = 1.14, the 
value of t/l is 0.37 anl the error is approximately 5 percent. In 
figure 3 the rolling-moment distribution is shown as a function of x/c. 
The moment distribution is zero for x/c < 0*90 and x/c > 1.25$ these 
values correspond to the downstream and upstream extremities of the 
aileron lea ding edge. Further, the moment distribution is discontinuous 
at the boundary of regions II and III. The chordwise pressure distri- 
butions for - J ■■■■■ = 0.2, 0.6, and 1.1 are shown in figure 4. 

2(c) 


Langley Aeronautical laboratory. 

National Advisory Committee for Aeronautics, 
Langley Field, Va., January 31, 1958- 
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APPENDIX A 
SYMBOLS 


c wing root chord 

distance from apex to position of maximum span of wing 


s 


x,y,z 


wing semi span 

Cartesian coordinates (x is in the stream direction) 
dummy variables of integration 


i = y + iz 

Z(x) spanwise ordinate of wing leading edge 


w -i 


Z _ ^(y) chordwise ordinate of wing leading edge 

f(x) spanwise ordinate of control-surface leading edge, 

f(x) = BLf (X - Xf) 


f’(x) 


_ df 


t(x) 


■y - dt 

dZ x =c 


spanwise ordinate of wing trailing edge 


( x i>yi) 

( x 2^2) 

u,v,w 


upstream coordinates of control-surface leading-edge comer 
downstream coordinates of control-surface leading-edge corner 
control- surface deflection 

nondimens ional disturbance velocities in x-, y-, 
and z-directions , respectively 


q> 


disturbance velocity potential 



z 
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- 

< 

ii 

'x' 

v« ■ 

y)ay 



* 

V(x;£) = v 

- iw 




L 

lift 




M x 

moment about x-axis 




°P 

pressure coefficient on the wing 




P(xj£) 

polynomial in f; 




A,B 

coefficients of polynomial 




<3oo 

stream dynamic pressure 





defined by equation (2.21) 




■ G(x;y) 

defined by equation (2.5) 



- 

g(ijUx)) 

defined by equation (2.11) 




H(|jl(x)) 

defined by equation (2.15) 




H(ftiUx)) 

defined by equation (5*8) 




S(x) 

evaluated f ran. equation (2.14) 




M(x) 

defined by equation ( 2 . 15 ) 




Q(x) 

evaluated frcm equation (5-6) 




N(x) 

defined by equation (5-7) 




K,E 

cangplete elliptic integrals of first and second 
respectively, with modulus k 

kind. 


K',E' 

complete elliptic integrals of first 
respectively, with modulus k' 

and second 

kind. 


F(e,k),E(e,k) incomplete elliptic integrals of 

respectively 

first and second kind, 


n(e,a 2 ,k) 

elliptic integral of third kind 



«r 

Z(0,k) 

Jacobian zeta function 
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A Q (0,k) Heuman's lambda function 



I imaginary part 

ft (xjy,Tj) defined by equations (B3) 

r(|}Ti) defined by equation (B9) 

2 K 

, _ n 

5 K' 

_ _ A F(0,k) 

5 K 

T' 



NACA TN \ 2\2 




q. = e‘^> 
q' = 

Subscripts : 

o values on z = +0 

x partial derivative with respect to x 

y partial derivative with respect to y 
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APPENDIX B 

EVALUATION OF G(x;y), g(|;l(x)), AND BELATED FUNCTIONS 


The solutions for the slender wing with symmetric spanwise loading 
involve elliptic integrals of the third kind in region II for c < c^_ 

and in region III. These Integrals which arise in the evaluation of 
the functions G(x;y) and g(|$2(x)) and in the expression for the 
pressure in region III are reduced to a form suitable for computation 
by their relation to theta functions. The functions fl(x;y,T]) and 
r(£$q), which are logarithms of theta functions, are readily evaluated 
in series form. These relations, as well as the reduction of the 
elliptic integrals which arise in the solution for symmetric spanwiBe 
loading, are given in reference 8, for example. 


The Functions Related to the Pressure 

With an interchange of the order of integration and with the use of 
partial fractions, the function G(x;y) given by equation (2.5) can be 
evaluated as 




<*n (bi) 


where F(0,k), E(0,k), and ll(0,a 2 ,k) are elliptic integrals of the 

first, second, and third kind, respectively, and 




a 2 = 


l 2 - t 2 
l 2 - q 2 


For y = t, 
complete. 


the argument 0 is 


£, and the elliptic integrals are 
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This elliptic integral of the third kind in equation (Bl) can he 
evaluated in terms of the Jacobian zeta function and a theta function. 

For 1 < a, 2 < 


n(0,a.2,k) - - | ~ ^ } g |?(0,k)Z(m,lO - O(xji))] 


where 


and 


CO 

0(x 5 y,Ti) = i log + V & Bin( 2 m a)sin( 2 mT) 

' ' “ 2 sin(cr - r) m sinh.(2mp) 


m=l 


1 + 


^(xjy^Ti) = | log 


i 2 -yg 

12 . .2 


1 - 


fi 2 - y 2 
l 2 . n 2 


(for x = c) 


(B2) 


(B3) 


sin cu = — = 
a 


,2 2 
* . - a 

z 2 _ t 2 


_ « F(cp,k) 
2 K 


T - * g(gU) 

2 K 


F o w 


£ K! 

2 K 


= e-2p 


q = e 


From the definition of the zeta function z(0,k) = E(0,k) - — F($,k) 

K 

and equations (B2), the combination of the functions G(xjy) and G(x;t) 
which arises in the expression for the pressure is 



k6. 
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0(l,y) -£&£l O^t) ~—r 

K *i(z 2 -t 2 ) 


0(icjy)_ O(l)t) - 0 


ij/^5 +z ^ ,!£) X l^ 2 ‘'> 2 Xi £ - t£ )’i e -' i+ z f t w o x (*;n)«(*;i)4') (r > t) 


(A) 


(y - 1 ) 


For x near c, the function fi(x;Tj) can he approximated hy Idle first 
term of equations (B3). 


The Functions Related to the Integral Equation 

With an interchange of the order of integration and with the use 
of partial fractions in equation (2.11), the integral on y can he 
expressed as the sum of two integrals, and the function g(i;l(x)) is 
evaluated as 




1 (x) 


T) 2 -t 2 


T Z 2 (x)-Z 2 Z 2 ( 


+^s(e,k') ♦^n(e,a'2,k') 


z 2 -t 2 


z 2 - n 2 


x)-t £ 


d?} 


(B5) 


where 


sin 8 = 


' Z 2 (x) - l 2 
l 2 (x) - t 2 


k- 

Z 


a 


' 2 = 1 - a 2 =-H 


2 _ t 2 


Z 2 - q 2 


From equations (Bl) and (B5), the quantity g(£;Z(x)) + G(g;t) 

K 

which arises in the integral equation for S(x) can he expressed as 
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8(8jl(x))+£k£l<j{i}t) -jrf- 


2 

■ 

I t(x)k 2 , 

Xl (l 2 -t a ) 

^[l 2 (x)-J 2 ][l 2 U)-t 2 ) 2 K 


l Vo^Ujn) |^[l 2 -t| 2 X T ! 2 -t 2 )n 4n+H(t;l(x)) (B6) 


where 


H(§;Z(x)) — | ^ v 0g (S;ti) a ^ 


,4 s - 1 2 

Z 2 - t 2 

h 2 - t 2 

Z2 - r[2 


? ' o — - t o g(e,a' 2 ,k') - F(e,k*) - 




dr} 


(BT) 


and Henman's lambda function (ref. 10) is defined by the relation 
A o (0,k) = §[e F(0,k') + K E(0,k') - K F(0,k')] 


The parameter a,' 2 is zero for T| » t; then, H(0,O,k') = F(0,k') and 
the function H(g;Z(x)) is zero. For t) > t, the elliptic integral 

H(0,a' 2 ,k') can be expressed. in terms of Henman's lambda function and 
a theta function as 


H(0,a' 2 k' ) = F(0,k' ) - I ||^?[F(0,k* )A 0 (a),k' ) - rUjuj] (b8) 


where 


r(6jTj) = F(0,k' ) + tan -1 <( 


00 

2 (-l) m+1 q. ,ni sin(2mr' )sinh {2m(p'-a')] 


m=l 


00 

1+2 ^ (-l) m q ,m cos(2mr' )cosh |2m.(p'-cr' )j 


L m=l 


(B9) 
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and 


, _ jt F(o> f k) r' - - F ( 6 ' k ') p' = - — q' = e -2 ?' 

0 ~ 2 K* 2 K* P 2 K' H 

Substituting from equations (B2) and (b 8) into equation (BT), making 
use of the definitions of Ao(<x>,k’) and Z(o>,k), and using Legendre’s 

relation EK' + E'K - KK' = ^ give the function H(|;Z(x)) as 


E(S,l(*» - f l «o 8 (Sil) [ F(e ' k K?K (m,k) - ( B1 °) 


For small values of 0, the functions F(0,k' ) and n(0,a,' 2 , k' ) are 


I 2/ 2 

approximated by sin 0 = \ 1 > x ^ and frcm equations (B7) and (B2), 

V Z 2 (x) - t 2 

the function H(£jZ(x)) for x near c becomes 




t£Hr* +zM 


dtj 


(Bll) 


An alternate form for x near c is obtained by making use of the 
appro xima tion n(e,a ,2 ,k') = 1 tan“ 1 (a tan 0) + 0(k' 2 ). The function 

C 6 

H(£;Z(x)) then becomes 


+ 




i jg^g t „-iP 2 (^)rj a , h 2 w-i 2 

I n2 


Z 2_ n 2 Vl 2( x ). t 2 


Z 2 (x)-Z 2 

Z 2 (x)-t 2 


Z(cD,k)l 


dti 


H(ft;l(x)) 
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APPENDIX C 

NUMERICAL SOLUTION OF THE INTEGRAL EQUATIONS 
FOR THE FUNCTIONS S(x) AND Q(x) 


The kernels of the integral equations for the functions S(x) 
and Q(x), which arise in the solution for the pressures in region II 
for c < are equivalent to those given hy Mirels for the flat-plate 

wing at angle of attack and for roll. The integral equations for the 
cambered wings differ from those for the flat plate only in the left- 
hand members . Since these equations can be inverted analytically in the 
neighborhood of x = c, the inversion in general can only be accomplished 
by numerical methods. The numerical method of solution presented in this 
appendix is based upon the method developed in reference 5 • 


Solution for S(x) 

Let the range of integration be divided into n intervals of 
length A and let the end points of the intervals be denoted 

by Xj_ = c + i A and the midpoints of the intervals by Xj_ = c + ^i - ^A 
where i = 0, 1, . . . n and x = x^ The integral equation for S(x) 
given by equation (2.1 1 !-) can then be expressed as 


M(x) = £ J 1 S(|)l* 


1=1 i -1 

With the notation 




k^I l(x) 


K 


^[l 2 (x) - Z 2 ][l 2 (x) - t 2 ] 


dl 


j 




- 21 A o (9 ' k) 7 . 
2 K 


= K x i) 



l (h) 


t 


i 




and with the use of the. mean-value theorem, the integral equation is 
approximated as 
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M ( x n)= Y, S ^) 
i=l 




Hie values of sfci), s(x £ ), . . . are found by successively 
setting n = 1, 2, . . . . The function is regular everywhere 

and vanishes for | = c and £ = and the integrals can be evaluated 
by Simpson's rule or other suitable methods. 


Solution for Q(x) 

Hie integral equation (5*6) can be Bolved by the same method as 
that given for the evaluation of S(x). Dividing the range of integra- 
tion c to x into n intervals gives equation (5-6) as 


N(x) 


£ J x± sum 


i=l i-l 




The values Q(xJ, Q(x 2 ) , . . . are obtained by successively setting 

n = 2 , • • • • 
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